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Energy Harvesting Transmitters that Heat Up: 
Throughput Maximization under Temperature 

Constraints 

Omur Ozel Sennur Ulukus Pulkit Grover 


Abstract 

Motivated by damage due to heating in sensor operation, we consider the throughput optimal 
offline data scheduling problem in an energy harvesting transmitter such that the resulting temperature 
increase remains below a critical level. We model the temperature dynamics of the transmitter as a 
linear system and determine the optimal transmit power policy under such temperature constraints as 
well as energy harvesting constraints over an AWGN channel. We first derive the structural properties 
of the solution for the general case with multiple energy arrivals. We show that the optimal power 
policy is piecewise monotone decreasing with possible jumps at the energy harvesting instants. We 
derive analytical expressions for the optimal solution in the single energy arrival case. We show that, 
in the single energy arrival case, the optimal power is monotone decreasing, the resulting temperature 
is monotone increasing, and both remain constant after the temperature hits the critical level. We then 
generalize the solution for the multiple energy arrival case. 

I. Introduction 

In many wireless sensor applieations, temperature increase caused by sensor operation has 
to be carefully managed. For example, wireless sensors implanted in the human body have to 
be designed such that the temperature due to their operation does not cause any threat for the 
metabolism. A line of medical research started by Pennes in 1948 [HI explores the temperature 
dynamics due to electromagnetic radiation in conjunction with heat losses to the environment 
and dissipation of heat in the tissue. In the context of sensors that communicate data, temperature 
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sensitivity varies depending on the type of tissue. For a given speeific tissue, it is reeommended 
that the temperature does not exceed a critical level, in order to prevent damage to the tissue. 
This necessitates careful scheduling of data transmission O. This problem arises in various 
types of body area sensor networks, see e.g., Q-jH and references therein. Finally, temperature 
increase in a sensor is a threat for the proper operation of the hardware itself (H-IIH. In this 
context, the electric power that feeds the amplifier circuitry has to be carefully scheduled so as 
to avoid permanent damage in the circuit. 

In order to obtain design principles with regard to temperature sensitivity of such systems, 
determining transmission schemes under a safe temperature threshold Tc is a useful objective. In 
this paper, we consider data transmission with energy harvesting sensors under such temperature 
constraints. Data transmission with energy harvesting transmitters has been the topic of recent 
research ifTOjl - lfTTll . In particular, throughput maximization under offline and online knowledge 
of the energy arrivals is considered in these references for single-user and multi-user energy 
harvesting communication systems. In lfT^ - lf22]| . this problem is investigated under imperfections 
such as battery energy leakage, charge/discharge inefficiency, and presence of processing costs. 

In the current paper, we aim to bridge physical heat dissipation with data transmission in energy 
harvesting communication systems. When the sole purpose is to maximize the throughput, the 
transmitter may generate excessive heat while utilizing the energy resource. In a temperature 
sensitive application, the heat accumulation caused by the transmission power policy has to be 
explicitly taken into account. In such a case, heat generated in the transmitter circuitry causes a 
form of “information-friction” |[2^ . We study the effect of this “friction” in a deadline constrained 
communication of an energy harvesting transmitter over an AWGN channel. For simplicity, we 
use transmit power as a proxy for hardware power. That is, we assume that the energy dissipated 
by the power amplifier dominates other energy sinks in the circuitry. More work is needed to 
understand full implications of communication circuitry’s energy in this context. Our formulation 
also relates to [|2^ in that the cumulative effect of heat generated in the hardware affects the 
communication performance. 

We determine the throughput optimal offline power scheduling policy under energy harvesting 
and temperature constraints. Our thermal model is based on a view of the transmitter’s circuitry as 
a linear heat system where transmit power is an input as in [[U, 0, 10 , IIMI . We impose that the 
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temperature does not exceed a critical level T^. Consequently, we obtain a convex optimization 
problem. We solve this problem using a Lagrangian framework and KKT optimality conditions. 
We first derive the structural properties of the solution for the general case of multiple energy 
arrivals. Then, we obtain closed form solutions under a single energy arrival. For the general 
case, we observe that the optimal power policy may make jumps at the energy arrival instants, 
generalizing the optimal policies in IfTOl . IfTTI . Between energy harvests, the optimal power is 
monotonically decreasing. We establish for the case of a single energy arrival that the optimal 
power policy monotonically decreases, corresponding temperature monotonically increases, and 
both remain constant when the critical temperature is reached. Then, we consider the case of 
multiple energy arrivals. We observe that the properties of the solution for the single energy 
arrival case are guaranteed to hold only in the last epoch of the multiple energy arrival case. In 
the remaining epochs, the temperature may not be monotone and the transmitter may need to 
cool down to create a temperature margin for the future, if the energy harvested in the future is 
large. We illustrate possible cases and obtain insights regarding the optimal temperature pattern 
in the multiple energy arrival case. 


II. The Model 

We consider an energy harvesting transmitter node placed in an environment as depicted in 
Fig. [U The node harvests energy to run its circuitry and wirelessly send data to a receiver. 

A. Channel Model 

The received signal Y, the input X, fading level h and noise Z are related as 

Y = VhX + Z (1) 

where Z is additive white Gaussian noise with zero-mean and unit-variance. In this paper, the 
channel is non-fading, i.e., h = 1. We use a continuous time model: A scheduling interval has 
a short duration with respect to the duration of transmission and we approximate it as [t,t + dt] 
where dt denotes infinitesimal time. In [t,t + dt], the transmitter decides a feasible transmit 
power level P{t) and | log (1 -f P{t)) dt bits are sent to the receiver, where the base of log is 
2. To be precise, the underlying physical signaling is in discrete time and the scalings in SNR 
and rate due to bandwidth and the base of the logarithm are inconsequential for the analysis. 
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Fig. 1. The model representing an energy harvesting wireless node placed in an environment that has constant temperature Te. 


B. Energy Harvesting Model 

As shown in Fig. |2l the initial energy available in the battery at time zero is Eq. Energy 
arrivals occur at times {si,S25---} in amounts {i?i, £'2, ■ ■ with sq = 0. We call the time 
interval between two consecutive energy arrivals an epoch. D is the deadline. Ei and Si are 
known offline and are not affected by the heat due to transmission. Let h{t) = max{k : Sk < t} 
and N be the number of energy arrivals in the interval [0, £) and by convention we let sjv+i = D. 
Power scheduling policy P{t) is subject to energy causality constraints as: 

nt 

/ P{r)dT < Y, Ei, yt e [0, D] (2) 

>^0 i=0 

C. Thermal Model 

In our thermal model, we use the transmit power as a measure of heat dissipated to the 
environment. In particular, we model the temperature dynamics of the system as follows: 

jT{t) = aPit) - h{T{t) - Te) + c (3) 

where P{t) is the transmit power policy and T{t) is the temperature at time t. Tg is the constant 
temperature of the environment that is not affected by the heating effect due to the transmit power 
level P{t). a and b are non-negative constants, c represents the cumulative effect of additional 
heat sources and sinks and it can take both positive and negative values. In the following, we 
consider the case of no extra heat source or sink, i.e., c = 0. 
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Fig. 2. Energy Ei becomes available for data transmission at time Si. D is the deadline. 

Our thermal model in ([3]) is intimately related to the thermal model in 0, lUl where hardware 
heating is modeled as a first order RC heat circuit. In particular, thermal dynamics of a power 
controlled transmitter due to its amplifier power consumption (see e.g., Il25l f could be modeled 
as in ([3]). We also refer the reader to Il24l for a related heating model. Our thermal model is 
also related to the well-known Pennes bioheat equation [[T]|. We assume, for simplicity, that the 
spatial variation in temperature is not significant and leave the general case of spatial temperature 
variations as future work. 

From dU), the solution of T(t) for any given P{t) with the initial condition T{t') at time t' is: 

T{t) = (aP(r) + bT^) dr + (4) 

By inserting f' = 0 in (jl]), we get (c.f. Il24l Eq. (3)]): 

T{t) = e-“ (aP(r) + hTe) dr + T(0)^ (5) 

The temperature should remain below a critical temperature Tc, i.e., T{t) < Tc, where we 
assume that Tc > Tg. Let us define Ts = — Tg, which is the largest allowed temperature 

deviation from the environment temperature. Typically, initial temperature is Tg, i.e., initially the 
temperature is stabilized at the constant environment temperature Tg. From ([5]), using T(t) < Tc 
and T(0) = Tg, we get the following equivalent condition for the temperature constraint: 

ae^^P{T)dT < Tse^\ Vt G [0, D] (6) 

Note that the temperature constraints in ® and the energy causality constraints in Q do not 
interact. Due to the heat generation dynamics governed by ([31), we observe in ® that the cost 
of power increases exponentially in time (i.e., the multiplier in front of T(t) is exponential in 
r) while the heat budget also increases exponentially in time (i.e., the upper bound on the right 
hand side of ® is exponential in t). 
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III. Problem Formulation 

Offline throughput maximization problem over the interval [0, D] under energy eausality and 
temperature eonstraints with initial temperature r(0) = is: 


r-D 


max 


p(i), telo,D] Jq 2 


log (1 + P{t)) dr 


s.t. 


/ ae^"F(r)dr < T^e^, Vf 

Jo 

/ P{T)dT<'^Ei, yt 
Jo ,_n 


(7) 


where the spaee of actions is the set of measurable functions P{t) defined over the interval 
[0, D]. Note that (|7]) is a convex functional optimization problem. 

The Lagrangian for (|7]) is: 


-D 


c = 


1 

2 

rD 


rD 


log (1 + P{t)) dt - I X{t) ( / ae^P{T)dT - Tse^^ ) dt 

m 


^ nt h{t) \ 

/ P{T)dT — 'S^Ei I dt 

/o ..0 J 


( 8 ) 


Taking the derivative of the Lagrangian with respect to P{t) and equating to zero: 


i-D 


i-D 


Jxt 


1 + P(t) 


\{T)dT 


j3{T)dT = 0 


(9) 


which gives 


p(t) = 


S° + e“ /y X{T)dT 


r-D 


( 10 ) 


In addition, the complementary slackness conditions are: 


\{t) / ae^^P{T)dT - Tse^^ =0, Vf 


/3(f) 


' Pt Mi) \ 


= 0, Vf 


( 11 ) 

( 12 ) 


In dH) and (fT]1)-(fT2l). A(f) > 0 and /3(f) > 0 are distributions that are allowed to have impulses and 


their total measure over [0, D] interval are not both zero, i.e., \{T)dT > 0 or l3{T)dT > 0, 
in order to prohibit P(f) from being unbounded. We note that (jll) and (fTT]) - (fT^ are necessary 


rD 
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and sufficient conditions since the problem is convex. The solution is unique almost everywhere 
as the objective function is strictly concave. 

We note that the problem in (|7]) could be solved by using calculus of variations. See |I71 
for application of calculus of variations for a similar problem to (|7]). As another alternative, 
we note that (|7]) could equivalently be solved by using a Hamiltonian approach from optimal 
control theory. In particular, we can cast the problem in ([7]) as an optimal control problem with 
pure state constraints lf26]l . In this case, the state of the system is the tuple [T(t) B{t)] where 
B{t) = fg P{T)dT is the total energy expenditure by the time t. The input is P{t) for Q <t < D. 
This problem is in the following form: 

rD ^ 

max / -\og(l + P(t)) dr 
Pit), mo, D] io 2 ^ 

s.t. jT(t) = MT, B, P), = MT, B, P) 

gi{T, B, t) < 0, g 2 iT, B,t) <0 (13) 

where fi{T,B,P) = aP - b{T - T^) and / 2 (T,S,P) = P while gi{T,B,t) = T - and 
g 2 {T,B,t) = B — Yli=o Note that gi and g 2 do not depend on the input P. With these 
selections, optimization problem (fT3l) is in the same form as that stated in ll2^ Eqs. (2.1)-(2.6)]. 
In this case, Hamiltonian is 

P(T, B, P, Ai, A 2 , t) = ^ log (1 + P) - Ai(f)/i(P, B, P) - A 2 (t)/ 2 (T, B, P) (14) 
and the corresponding Lagrangian is 

CHiT,B,P,Xi,X2,t) = 'H{T,B,P,Xi,X2,t) - ui{t)gi{T, B,t) - U2{t)g2{T, B,t) (15) 

where Ai(t) and X 2 {t) are the co-state trajectories; ui{t) and mit) are multiplier functions. We 
note that Pontryagin’s maximum principle is necessary and sufficient in this case since (flSl) is 
a concave maximization problem. One can derive the equivalence of necessary and sufficient 
conditions for this optimal control problem to those in @ and (fTT]) - (fT2l) . 

In the following, we proceed with the Lagrangian formulation in ([8]) and the corresponding 
optimality conditions in (|9l) and (fTTI)-(fT2l). 


IV. General Properties of an Optimal Policy 


In this section, we obtain the structural properties of the optimal power scheduling policy 
using the optimality conditions. In the following lemmas, P{t) refers to the optimal policy and 
T{t) is the resulting temperature unless otherwise stated. 

We first note that the temperature level never drops below Tg. In particular, if the initial 
temperature is between Tg and Tg, the temperature at all times will remain between Tg and Tg. 

Lemma \ Tg < T{t) < Tg whenever the initial temperature is Tg < T{0) < Tg. 

Proof: From dH), since P{t) > 0 we have ^T{t) > 0 whenever T{t) = Tg. The constraint 
T{t) < Tg is satisfied by any feasible policy in (|7]). ■ 

The following lemma states that if the temperature T{t) is constant, then the power P{t) is 
constant also (while it is not true the other way around, see Lemma [3]), and that if the temperature 
hits the maximum allowed level Tg, then the power must be below a threshold. 

Lemma 2 Whenever T{t) is constant over an interval I C [0, D], P(t) is also constant over 
that interval. If the temperature hits the level Tg at t = th, then P{th + e) < ^ for all sufficiently 
small e > 0. 

Proof: If T(t) is constant in I, jfT{t) = 0 and from dS]), P{t) is also constant in the same 
interval. If Tfh) = Tg for some th G [0, D), then ^T{th + e) < 0 and from dl]), P{th + e) < 

■ 

The following lemma shows that if the power P{t) is a monotone increasing function, then 
so is the temperature T{t). We first prove this result for piecewise constant functions and then 
generalize it to arbitrary functions. We note that a particular instance of a monotone increasing 
piecewise constant power is observed in the solution of the throughput maximization problem 
without temperature constraints ifTOl . 

Lemma 3 If P{f) is a monotone increasing piecewise constant function, then T{t) is monotone 
increasing. More generally, if P{t) is a monotone increasing function, so is T{t). 

Proof: We first prove the first statement of the lemma which is concerned with piecewise constant 
functions. Let us start with the case of a single constant power value for the entire duration of 


communication, i.e., P{t) = p for t e [0, D]. From Q, we have: 


T{t) = ( f {ap + bTe) dr + T(0) 


( 16 ) 


= e 


Te + ^P+ (t( 0) - T, - ^p) e-“ 


(17) 


(18) 


For T(0) = Te, (fTSl) is a monotone inereasing funetion of t. In partieular, T{t) < Tg + |p. Now, 
let us consider the ease of M constant power levels for the duration of eommunication, i.e., 
P{t) = Pi over the interval [/j-i, h) where Pi < p^+i for all i and 0 = lo < h < ■ ■ ■ < Im — D 
where M > 1 is the number of intervals. In this ease, we have for t G [li-i, h): 



(19) 


where T(/j_i) < Tg + fpi-i- Hence, the coefficient of in (fT9l) has a negative sign as 

T(Jj_i) — Tf. — ^pi < j {pi-i — Pi) < 0. This proves that T{t) is monotone inereasing. 

To generalize this result for any monotone inereasing funetion P{t), we obtain any monotone 
inereasing simple approximation [l27ll of P{t), denoted as Pn{t), such that Pi{t) < P 2 {t) < 
... < Pnit) for all t G [0,79] and Pnit) P{t) pointwise. For example, one ean select 

Pn{t) = P(4(i_i)) for t G [Inii-i),Ini) and ^ (z - 1) for i = 1, ..., 2”. Let us call the 

resulting temperature Tn{t). Henee, e^*Pi(f) < e“P 2 (f) < ... < e'^^Pnif) for all t G [0, P] and 
e^*P„(f) e^^P{t) pointwise. By monotone eonvergenee theorem IfTTlI . we have 



( 20 ) 


Accordingly, r„(f) —)■ T{t) pointwise and we have 


jTn{t) = aPn{t) - b {Tn{t) - Pg) ^ jT{t) = aP{t) - b {T{t) - Pg), Vf G [0, D] (21) 


Sinee Pn(f) is a monotone inereasing pieeewise eonstant funetion, from the first part of the proof, 
Tn{t) is monotone increasing, i.e., ^Pn(f) = aPnit) — b{Tn{t) — Pg) > 0. Sinee ^Pn(f) ^ 
^P(t) pointwise, this implies ■^T{t) > 0, i.e., T{t) is monotone inereasing as well. ■ 

The next lemma shows that if the temperature remains eonstant over an interval, then that 


level eould only be Pg or Pg, i.e., any other temperature eannot be a stable temperature. 
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Lemma 4 If T{t) is constant over an interval I C [0, D\ then that constant level could only 
be Tg or T^. 

Proof: Assume T{t) is constant over I. Without loss of generality, assume that there is no energy 
arrival in the interval I, and otherwise let I be the portion of the interval without any energy 
arrivals. By Lemma [2l P{t) is constant over I. If P{t) = 0 over I, then T{t) = Tg from ([3]). If 
P{t) f 0, we have from (fT^ 

P{t) = —n -^-n-1 (22) 

/?(r)(ir + X{T)dT 

where (3{t) = 0 over the interval / by (fT^ since f3{t) >0 implies energy constraint is tight and 
P{t) = 0. Therefore, (3 [t) dr = S is constant over 1. If T{t) < Tc, then by (fTTl) . X{t) = 0 
over I and hence X{T)dT = C is constant over I. However, this makes (l22l) a time varying 
function of t because of the term in the denominator, and this contradicts the fact that P{t) is 
constant. Finally, if C = 0, this means that the temperature constraint is never tight. In this case, 
the piecewise constant power policy in ifTOll is optimum, and the temperature is monotonically 
increasing from Lemma [H and therefore, cannot be a constant over an interval. ■ 

The following lemma states that at the end of the communication session either the harvested 
energy is exhausted or the critical temperature is reached. 

Lemma 5 At t = D, either the temperature constraint or the energy causality constraint or 
both are tight. 

Proof: If neither of the constraints are tight, then the power policy P{t) could be increased over 
a set of non-zero Lebesgue measure in the last epoch. This strictly increases the throughput, 
contradicting the optimality. ■ 

The following lemma shows that the optimal power should be monotonically decreasing 
between energy harvests. 

Lemma 6 P{t) is piecewise monotone decreasing except possibly at the energy arrival instants. 
In particular, it is monotone decreasing between consecutive energy harvests. 

Proof: We prove the statement by contradiction. Assume that for some interval P{t) 

is strictly monotone increasing, and that the interval [^ 1 ,^ 2 ] does not contain an energy arrival 



11 


P{'r)dr 

instant. Define a new power policy as Pnewif) = ^ ^ [^ 11 ^ 2 ] and Pnewit) = P{t) 

otherwise. Pnew{t) satisfies the energy causality constraint in (|7]) since Pnewit) uses the same 
amount of energy as P(f) over [ti,t 2 ] and the energy constraint for Pit) is not tight in this 
interval. Pnewit) also satisfies the temperature constraint. To see this, we first note that Pnewit) 
satisfies the following inequality (see ll28l Theorem on p. 207]): 


rt2 


ae’^^Pn 


rt2 


}i'r)dT < / ae°'^PiT)dT 






(23) 


as both Pit) and are monotone increasing. In addition, since P(f) is temperature feasible: 



ae^^PiT)dT < 


t2 


ae^"P(r)dr < 


(24) 

(25) 


Jo 

Combining (l23l) and (1251) . we conclude that Pnewit) satisfies the temperature constraint at 1 = ^ 2 : 


^2 


I't2 


ae^'" Pnewir)dT = / ae^'^ Pnewir)dr + / ae^^P„e«,(r)(ir 


rti 


< I ae^^P(r)(ir + / ae^'^P(r)(ir 


rt2 


br ] 


fti 


< Tse 


bt2 


(26) 

(27) 

(28) 


Additionally, the temperature constraint is satisfied for t > t 2 since Pnewit) and P(f) are identical 
for t > t 2 and P(f) is temperature feasible. Hence, we need to show that Pnewit) satisfies the 
temperature constraint for all t G to establish the temperature feasibility of Pnewit)- That 

is, we need to show: 


ae^'"P(r)dr+ [ Pnewi'T)dT <Tse^\ fe(fi,f 2 ) 


'0 


Since Pnewit) = p is constant over [ti,t 2 ], we have: 

/ ae'"'^Pnewir)dr = % (e“ - , t e [ 11 ,^ 2 ] 


Using (1301) in (1291) and since e’^ > 0, (1291) takes the following equivalent form: 


^—bt 


ae'’'^P(r)dr - ] +^P<Ts 


a \ a 


(29) 


(30) 


(31) 
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Note that the left hand side of (l3T]) is either monotone inereasing or monotone decreasing in 
t as it is a linear function of e~^^. Since the inequality (|3TI) holds at t = ti and t = t 2 as 
Pnew{t) satisfies the temperature constraint at those points, we conclude that Pnew{t) satisfies 
the temperature constraint for all t E [^ 1 ,^ 2 ]- In addition, Pnew{t) yields higher throughput than 
P{t) due to the concavity of logarithm. This contradicts the optimality of P{t). The proof holds 
even when [ti, ^ 2 ] includes an energy arrival instant provided that the energy causality constraint 
is not tight at that instant. ■ 

Next, we show that discontinuities in the power level could only occur in the form of positive 
jumps, and only at the instances of energy harvests. 

Lemma 7 If there is a discontinuity in P{t), it is a positive jump and it occurs only at the 
energy arrival instants. The temperature T{t) is continuous throughout the [0, D] interval. 

Proof: Since is a continuous function of t, X{t) > 0 and (3{t) > 0, any jump in P{t) has to 
be positive due to (fTOl) . Any positive jump at instants other than Sk violates monotonicity of P{t) 
within each epoch due to Lemma 0 Due to ([5]), the resulting temperature T{t) is continuous 
throughout the [0, D] interval. ■ 

By Lemma|71 we can take /3(f) in the form (3(f) = (3j6(f—Sj) without loss of optimality, 

where I3j > 0, j = 1,..., N + 1, are finitely many Lagrange multipliers corresponding to the 
energy causality constraints at the energy harvesting instants sj and the deadline, sn+i = D. 

The next lemma shows, for an arbitrary feasible policy P{f), that if the temperature reaches 
the critical level Tc at some th, then the power just before th must be larger than a threshold. 

Lemma 8 If Tiff) = Tc for some th E [0, D), then Pifh — e) > ^ for all sufficiently small 
e > 0 . 

Proof: Since Tifh) = Tc, we have: 



(32) 


We combine ® with (l32l) to get 



Vf e [ 0 , 4 ] 


(33) 
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which implies in view of the eontinuity of P{t) (exeept for the finitely many energy arrival 
instants) proved in Lemma |7] that P{th — e) > ^ for all sufficiently small e > 0. ■ 

We next state the eontinuity of the optimal power poliey P{t) at points when it hits the eritieal 
temperature Tc. 

Lemma 9 IfT{th) = Tc for some th G [0, D) then P{t) is continuous at th and Pith) = 

Proof: The proof follows from Lemma [2] and Lemma [8] and the fact that negative jumps in P(f) 
are not allowed due to Lemma IT] ■ 

Next, we show that when the temperature hits the boundary it has to return to Tc. 

Lemma 10 Whenever = Tc for some th < D, there exists t > th such that Tit) = Tc. 

Proof: Assume that Tfh) = Tc for some th < D and T(f) < Tc for all th < t < D. 'Qy 
Lemma m Pith) = From (jH) with Tith) = Tc, the eonstraint T(f) < Tc beeomes: 

[ ae^"F(r)dr < Ts (e“ - , th < t < D (34) 

Jth 

Sinee T(f) < Tc in th < t < D, only energy eausality eonstraint is aetive and thus P(f) for 
th < t < D is the pieeewise eonstant monotone power alloeation in ifTOll . On the other hand, 
Pit) = ^ satisfies (l34l) with equality for all t. Therefore, we must have P(t) = c < ^ for all 
t ^ ith,th + <^) for some () > 0. However, this eontradiets Pith) = ^ sinee there eannot be a 
negative jump in P(f) by Lemma |71 ■ 

The following lemma identifies the exact conditions where the power P(f) makes a jump. 

Lemma 11 If there is a jump in Pit), it occurs only at an energy arrival instant, when the 
battery is empty and the temperature is strictly below Tc. 

Proof: Due to the slaekness eonditions in (fTT]) - (fT^ . a jump oeeurs if either the battery is empty 
or the temperature eonstraint is tight, i.e., P(f) = Tc. By Lemma HI P(f) is eontinuous whenever 
T(f) = Tc. Therefore, a jump in P(f) oeeurs at an energy arrival instant, when the battery is 
empty and P(f) < Tc. ■ 

We finally remark that energy may have to be wasted as aggressive use of energy may eause 
temperature to rise above the critical level. 
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V. Optimal Policy in the Single Energy Arrival Case 


In this section, we consider a single epoch where E units of energy is available at the 
transmitter at the beginning. We first develop further structural properties for the optimal power 
control policy in this specific case and then obtain the solution. 

A. Properties of an Optimal Policy 

The next lemma shows that, if the power falls below a certain threshold at an intermediate point 
and remains under that threshold until the deadline, then it should remain constant throughout. 

Lemma 12 T/'O < P{f) < ^ for t G [ti, D], then P(t) is constant over [fi, D\. 

Proof: Assume P{f) is not constant over [ti,D]. Let Er = P{t) dr > 0. Define a new 
policy Pnewif) = for t G [ti,D] and Pnew{t) = P{t) otherwise. Pnew(t) is both energy 
and temperature feasible. Energy feasibility holds by construction as Pnew and P have the same 
energy over [ti,D]. Temperature feasibility also holds: Tfi) < Tc since P{t) is temperature 
feasible and as we have T{t) < Tc for all p < t < D from ®. Now, by 

Jensen’s inequality Pnew{t) achieves strictly larger throughput since log is strictly concave. This 
contradicts the optimality of P{f). Hence, P(f) = c > 0 for f G [p, D].M 

The following lemma states that the power has to remain constant at the level ^ when the 
temperature reaches the critical level Tc. 

Lemma 13 Let t' G [0, D] denote min{f G [0, D] : T{t) = Tc\. If t' < D, then P{t) = ^ for 
all t G [f, D]. 

Proof: By Lemma [9l P{d) = ^. By Lemma P{t) is monotone decreasing, and thus 0 < 
P{t) < ^ for t' < t < D. By Lemma \T2\. P{t) = c for all t G [P, D]. By Lemma |71 P(t) is 
continuous and therefore, P{t) = ^ for all t G [P, D]. ■ 

The following lemma states that the optimal power is always larger than a constant value 
determined by the fixed system parameters. 

Lemma 14 The optimal policy P{t) satisfies: 



(35) 
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Proof: If the temperature constraint is not tight, then the problem reduces to the energy con¬ 
strained problem in which case P{t) = ^. If the temperature constraint is tight, P{t) is monotone 
decreasing by Lemma [6] and when the temperature level reaches Tc, P{t) remains at ^ by 
Lemma [T3l Hence, P(t) ■ 

The following lemma shows that, since the power is always larger than a constant value, 
battery energy level is never zero, except possibly at the deadline. 

Lemma 15 In an optimal policy, energy in the battery is non-zero except possibly at t = D. 

Proof: By Lemma [T4l the optimal power is always larger than a positive constant. Thus, the 
battery energy does not drop to zero. ■ 

The following lemma shows that the temperature is monotone increasing throughout the 
transmission duration, and also is a concave function of time. 

Lemma 16 The temperature with the optimal power policy is monotone increasing and concave. 

Proof: If the temperature constraint is never tight, then the optimal power level is and from 
Lemma [3l the temperature is monotone increasing. Concavity in this case follows from the 
concavity of the explicit expression in (fT^ with r(0) = Tg. Now, assume that the temperature 
constraint is tight at t = D. By Lemma [T4l P{t) > From ([3]), we have: 


= aP{t)-h{T{t)-Tc) 

(36) 

>J^-b{T{t)-Tc) 

a 

(37) 

= b{Tc-T{t))>0 

(38) 


as T{t) < Tc by the temperature constraint. Since P{t) is monotone decreasing by Lemma [6] 
and T{t) is monotone increasing, from (l3^ . ^ is monotone decreasing, proving the concavity 
of T{t) in this case. ■ 

B. Optimal Policy 

In view of Lemma [l5l the energy constraint can be tight only at t = D. Therefore, the 
corresponding Lagrange multiplier is a single variable (3{t) = — D). From Lemma [HI T{t) 
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is monotone increasing. Due to Lemma [T3l when T(t) reaches Tc, power level has to remain at 
Accordingly, we denote the instant when the temperature reaches Tc as to. 

1) Sufficiently Large Energy: In this case, the energy constraint is never tight, and /3 = 0. In 
view of Lemma [51 the temperature constraint is tight dX t = D. 

First, consider the case that D is sufficiently large so that there exists to < D such that 
T{to) = Tc- For t e [0,fo), T{t) < Tc and from (fTTl) . \{t) = 0. From (ITOl) . when t e [0,fo) we 
have P{t) = — 1 where C = X{T)dT > 0. Since at t = to the temperature reaches Tc, 

from Lemma [T3l we have P{t) = ^ for f e [to, D]. Then, the optimal power has the form: 

P(t) = - 1 


{uiO - uit - to)) + —u(t - to) 

a 


(39) 


where u{t) is the unit step function. Now, from Lemma |9l P{t) is continuous at to and C should 


be chosen accordingly. In particular, C = 


verifies 




-Wo following Lagrange multiplier \{t) 


b 

^ I + ^Tsb I ~ 

The corresponding optimal temperature pattern for 0 < t < fo is: 


(40) 


T{t) = a 



te-^d-to) + « -w - ^ + Te 

b b 


(41) 


and T(t) = Tc for to < t < D. We note that to satisfies: 







(42) 


so that T{to) = Tc- Hence, T{t) monotonically increases till it reaches which is consistent 
with Lemma [T6l 

Next, consider the case that D < to- In this case. 


Pit) = -e-“ - 1 


(43) 


where C = 


D 




Yy and A(f) = CSit — D). Therefore, the optimal P(f) in this case is 




^<5 , i-A .61) 3 \ -bt 


5'" 


(44) 
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We also remark that to level that satisfies (l42l) monotonieally increases with T^. To see this, 
we rearrange (1^ as follows: 


1 

b 




— to 


0 


(45) 


Let us define a multi-variable real function w{to, Ts) as the left hand side of (1451) and denote a 
specific solution as tg for fixed Ts. It is easy to see that (l45l) always has a solution to for fixed 
Ts. To see this, we evaluate the derivative with respect to to as: 


^yj(to, Ts) = “o - 1 < 0, Vto > 0 (46) 

dto + 1 ) 

That is, w{to,Ts) is monotone decreasing with to- At to = 0, w{to,Ts) > 0 while w{to,Ts) 
—oo as to grows. In view of the continuity of iv{to, Ts), there exists a to such that w{to, Ts) = 0. 
Additionally, we observe in (l45l) that for fixed to, w{to,Ts) monotonieally increases with Ts. 
Therefore, if w{tQ,Ts) = 0, then, due to monotone increasing property with respect to Ts, 
w{tQ,Tg) > 0 for Tg > Ts. Hence, for tg* such that w{t*Q ,T'^) = 0, we have tg* > tg due to 
monotone decreasing property with respect to to- 

2) Energy Limited Case: Note that the optimal power policies in the energy unconstrained 
cases in (l39l) and (l44l) have finite energies. If the available energy E is larger than the corre¬ 
sponding energy level in (1^ and (l44l) . then the solution is as in (1^ and (l44l) . Otherwise, the 
energy constraint is active and the Lagrange multiplier is /3 > 0. From (fT^ . we have: 


m = 


/3 4 - X{T)dT 

We first note that there is a critical energy level Ecritical such that if i? < Ecriucah then 


- 1 


(47) 


constant power policy P{t) = ^ is optimal. This critical level is: 


^critical 


Tsb De 


bD 


a e 


bD 


(48) 


This is the critical level below which the temperature constraint is not tight by the constant power 
allocation P{t) = The expression in (l48l) is evaluated from (fTSl) by inserting T(0) = Tg, and 
requiring T(D) < Tc. When E < Ecruicai, X{t) = 0 since temperature constraint is never tight. 
In this case, (3 = -e^- Ecritical is the maximum energy level for which a constant power level 
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is optimal. If P{t) = x{t) is monotone increasing over [0, D] and reaches Tc at f = D. 

If E > Ecriticah the constant power level does not satisfy the temperature constraint. We 

note from (|4^ that Ecritical increases with the deadline D. Therefore, there exists a deadline 
level D for which D > D implies E < Ecruicai and hence constant power policy is optimal. 

An alternative way of observing the behavior of the optimal policy is to fix the available 
energy E and and vary the critical temperature Tc. In this case, there is a critical temperature 
limit for which P{t) = ^ is optimal whenever 


_ rjn I ^ 


a E — 1 


" ' hD 


(49) 


which again is evaluated from (fTSl) with T(0) = Tg. In the following, we consider E > E, 


critical 


or Tc < so that both energy and temperature constraints are tight at the end of the 

communication session. 

Again, we consider two possibilities: temperature constraint becomes tight at a to < D, and 
temperature constraint becomes tight at t = D. In both cases, the energy constraint becomes 
tight at t = D. 

First, consider the case that to < D: Due to (fTTl) . X{t) = 0 for f e [0,fo) and from (fTO . we 
get: 

1 


P{t) = 


/3 + 


- 1 


(50) 


where C = \{T)dT > 0. Additionally, P{t) = ^ for the remaining portion of the epoch in 
view of Lemma [T3l to is such that for t > to, P{t) = ^ and T{to) = Tc. Since P{to) = ^ 
we have: 


T,b 


+ 1 


/3 + a 

Similarly, for T{to) = Tc, we have from ([5]) with r(0) = T^. 

1 


(51) 


1 I + bTc \ dt + Tc\ — Tf 


Finally, the energy constraint has to be satisfied at t = D\ 




(52) 


(53) 
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P(t) 


Pit) 




(b) Fixed E and varying Tc. 


Fig. 3. The optimal power policy in the single energy arrival case for different energy and deadline constraints. 


If there exists < D for (|5TI)-(|5^. then P{t) is: 

In this ease, the eorresponding Lagrange multiplier is: 

\{t) = bCe-^^^-^°^u{t - to) + - D) (54) 

Otherwise, when no such to < D exists, the temperature constraint is tight only at to = 

In this case, P{t) is as in (l50l) for f e [0, D] where jS and C have to satisfy: 

The corresponding Lagrange multiplier is \{t) = CS{t — D). 


Depending on the energy E and the critical temperature Tc, the optimal power scheduling 
policy P{t) varies according to the plots in Fig. [3l For small E and fixed Tc or for large Tc 
and fixed E, a constant power policy is optimal. For moderate and large E, the optimal power 
policy is exponentially decreasing and may hit the power level Note that to level at which 
temperature touches the critical level decreases as Tc is decreased and as E is increased. In 
particular, for fixed Tc, the level of to is bounded below by the solution for i? = cxo whereas for 
fixed E, to goes to 0 as Tc approaches Tg. 
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VI. Optimal Policy for Multiple Energy Arrivals 

In this section, we extend the solution to the case of multiple energy arrivals. We start with 
extending the properties observed for the single energy arrival case when initial temperature 
T(0) is different from T^,. The following lemma generalizes Lemmas [T3] and [T^ for the case 
of an arbitrary T(0). 

Lemma 17 Assume that the initial temperature T(0) is in the range Tg < T(0) < Tg instead 
o/r(0) = Tg and consider the single energy arrival case: P(t) is monotone decreasing. Let 
th e [0, D] denote min{t G [0, D] : T{t) = Tg}. Ifth < D, then P{t) = for all t G [th, D] and 
the temperature is monotone increasing and concave. IfT{0) = Tc, then P{t) = min hi¬ 
proof: If T(0) is in the range Tg < T(0) < Tg then, instead of db]), we have the following 
temperature constraint: 

ae^^P{T)dT < T^e^ - Tg, Vf G [0, D] (57) 

where Tg = T(0) — Tg > 0. Note that T^e^ — T^ > 0 for all t G [0, D], i.e., the right hand 
side of (1571) is always non-negative. The argument in Lemma is valid in the presence of the 
additional term Tg in (l57l) . and therefore P{t) is monotone decreasing. 

The second claim follows from the argument in Lemma [l3l In particular, in addition to 
Lemma IS Lemma \V2\ directly extends with the constraint in (l57l) . Hence, the result follows by 
applying the argument in Lemma [T3l 

Linally, T{t) is monotone increasing and concave due to the steps followed in Lemma ITS In 
particular, if the temperature constraint is tight at f = T, P(t) > Hence, (l3S) - (l38l) hold and 
the temperature is monotone increasing and concave. If T(0) = Tg, then P(t) = min 
due to the energy constraint. Note that the temperature decreases in case min = %■ ® 

As in the single epoch case, we will investigate the solution under special cases. In particular, 
we will investigate the solution according to the time when the temperature hits the critical 
level. To this end, we specialize in an interval [ti, ^ 2 ] such that T{t) < Tc for all f G [fi, ^ 2 ) and 
T(f 2 ) = Tc where 0 < p < t 2 < D. Note that the temperature T(f) is a continuous function of 
t and hence there exist such intervals. We assume that the solution is known in [0, p) U (^ 2 , D] 
and we let Tc < T(ti) < Tc. In this case, the solution of dV]) over the interval [^ 1 ,^ 2 ] is equal 
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to the solution of the following problem obtained by restricting the temperature constraint to be 
satisfied at f = ^2 only: 


max 

p{t), te[ti,t2] 

s.t. 


^ log (1 + P(t)) dr 




ae^^P{T)dT = - Tg 

h[t) 


f P{T)dT<^%, 

■^0 i=0 


yt E [^ 1 ,^ 2 ] 


(58) 


where Tg = T{ti) — > 0. In (1581) . Ei is determined as follows: Eq is the available energy in 

the battery at time t = p. Ei for i = 1,..., iV are the energy arrivals at instants Sj E (^ 1 ,^ 2 )- 
h(t) is defined accordingly. While the times Si are exactly those in the original problem, the 
amounts Ei may be different from the original amounts as some energy may be left for use in 
the {t 2 ,D] interval. For the following argument, whether Ei equals the original energy arrival 
amount is not relevant and we leave Ei as arbitrary amounts. To obtain the solution of (1581) 
using this Lagrangian framework, it is necessary and sufficient to find N + 2 variables jdi > 0, 
^ = 1,..., iV + 1 and C > 0 such that 


p(() = 


LES‘ft+Ce" 


t E [si-1, Si), z — 1,..., iV + 1 


(59) 


with the corresponding slackness conditions. Therefore, for the [ti,t 2 ] interval, the solution has 
the structure in (I59l) . which is parameterized by finitely many Lagrange multipliers. In particular. 


throughout an epoch over which T{t) < T^, power level satisfies P{t) = 




- 1 


for some 


j3 > 0 and (7 > 0 not both equal to zero. This also holds in a subinterval of an epoch over 
which T{t) < Tc. In the following lemma, we show that in such an epoch, the temperature T{t) 
is unimodal. 


Lemma 18 If P{t) = 


p+Ce^ 


- 1 


T{t) is unimodal over \ti,t 2 \. 


for t E [^ 1 ,^ 2 ] for some (3 > D and (7 > 0, the resulting 


Proof: From dH), we have for t E [p, ^ 2 ], 


T(t) = e 






_I3 + 


- 1 


+ bTf, ] dr + T{ti) 


(60) 
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First, we note that when P{t) = 0, -^T(t) < 0 from Q- Hence, it suffices to show that T{t) is 
unimodal when P{t) = — 1 > 0. By evaluating the integral, we get 


a 




—bt 




(61) 


We claim that T{t) in (IMl) is unimodal for t > p. Note that the derivative of T{t) is: 




f ae^* a 

V/S + C'e“ “ C 


(13 + \ 

V/9 + Ce“i ) 


- h {Tit,) -T, + l) 



(62) 


We let X = x, = and concentrate on ^ log for x > x,. We note 

that ^ log ^ ^+cx ^ ) ^ strictly monotone decreasing function of a; for a: > Xi > 0. In 

particular, we have: 


d f ax (I3 + Cx\\ —Cx 

dx\(3 + Cx C^^\(3 + Cxi}} if^j^ Cxf 

Thus, log is strictly monotone decreasing in t. As log > 

0 at t = fi, we conclude that the factor in (l62l) that multiplies can take value 0 at most once. 
In particular, - § log - b (Tit,) - Te + f) can take positive or negative 

values at f = ti. If it is positive at f = fi, it hits value 0 at most once for f > fi. If it is negative 
at f = ti, it stays negative throughout t > t,. This proves that T(t) is unimodal over ® 

In the following lemma, we show that, in an epoch [sj,Sj+i], the temperature cannot return 
to Tc if it hits and falls below T^. 


Lemma 19 If Tith) = Tc and Tith + A) < Tc for some A > 0 where both th and th + A are 
in [sj, Sj+i], then Tit) < Tc for all t E [th + A, Sj+i]. 

Proof: By Lemma |9l Pith) = ^. By Lemma [6l power is monotone decreasing in an epoch. 
Therefore, if r(f/j+A) < Tc, then Pith+A) < ^ and hence Pit) < ^ for all t G [th+A, s^+i]. 
This, in turn, means that Tf) < Tc for all t E[th + A, s^+i]. ■ 

Next, we complete the unimodal structure of the temperature by showing that it has to be 
monotone decreasing if it hits and falls below Tc. 

Lemma 20 In an epoch [sj,Sj+i], if the temperature touches Tc at th and falls below it, then 
the temperature is monotone decreasing in [f/j,Sj+i]. 
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Proof: By Lemma [T9l if T(4 + A) < T^, then T{t) < Tc for all t ^ [th + A, Sj]. Therefore, we 
have 


Pit) 


1 

/3 + Ce“ 



t E [th + A, Sj+i] 


(64) 


for some (5 > 0 and (7 > 0. By Lemma [T^ T{t) is unimodal over t E [f/i + A, Sj]. Therefore, 
T(t) is monotone decreasing. ■ 

We next consider epochs [sj, Sj+i] and its subintervals over which T{t) < Tc and T{t) = Tc. 
By Lemma [T^ and in view of the discussion around (1581) . whenever T{t) < Tc over an epoch, 
T{t) reaches its peak level over that epoch at only one instance. Consequently, if T{t) < Tc 
for all t G [sj,Sj+i], there are three possible cases. The first two possibilities are that T{t) is 
monotone increasing or monotone decreasing throughout the epoch. The third possible case is 
that T(t) is monotone increasing in and monotone decreasing in for some 

tii G (sj,Sj+i). Otherwise, T{t) hits Tc and T{t) does not return to Tc if it falls below it due 
to Lemma [T9l Therefore, if T(t) hits Tc in an epoch [sj,Sj+i], then that epoch is divided into 
three successive subintervals with la = Ii 2 = [tii,ti 2 ) and 1*3 = [ti 2 ,Si+i] for 

some Si < til < ti 2 < Si+i. T{t) is monotone increasing over la, remains at Tc over 1*2 and is 
monotone decreasing over 1^3. We finally note that if T(t) < Tc at t = D, then T{t) < Tc for all 
f G [0,i7]. This follows from Lemma [TO In this case, the temperature constraint is never tight 
and the optimal power policy is identical to the one in [fTOll . 

In Fig. m we plot the optimal energy expenditure for different values of critical temperature 
level Tc. We observe that as Tc is decreased, the temperature budget shrinks and the temperature 
constraint becomes more likely to be tight. In this case, energy is spent faster not to create 
high amounts of heat in the system. In general, there is a tension between causing unnecessary 
heat in the system and maximizing the throughput. While we have fully characterized this 
tension in the single energy arrival case, it needs to be further explored in the multiple energy 
arrivals case. In particular, when a high amount of energy arrives into the system during the 
progression of communication, the transmitter has to accommodate it by cooling down and 
creating a temperature margin for future use. While maximizing the throughput generally requires 
using the energy in the system to the fullest extent, the transmitter may have to waste energy due 
to the temperature limit. We investigate this tension in numerical examples in the next section. 
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Fig. 4. Energy expenditure with the optimal power policy with multiple energy arrivals. In view of the temperature constraint, 
as Tc is decreased, the energy is spent faster subject to energy causality. 

VII. Numerical Results 

In this section, we provide numerical examples to illustrate the optimal power policy and the 
resulting temperature profile. For plots in Figs. [51 [6l |7] and [H we set a = 0.1, b = 0.3, = 37 

and Tc = 38. Therefore, the critical power level is ^ = 3. 

In Figs. [5]and[6l we consider the energy unlimited scenario. In this case, the solution of (l42l) 
is found as fg “ 2.993. In Fig. [51 we set D = 2 < to and we observe that the optimal power 
policy is always above the level In this case, power strictly monotonically decreases while 
temperature strictly monotonically increases with temperature touching the critical level Tc at 
the deadline. In Fig. [^ we set the deadline as D = 3.5 > to. We calculate that the energy 
needed to have the power policy in Fig. [^ is = 17.98. In other words, if the initial energy 
is > 17.98 then the power policy in Fig. [^ is optimal. We observe that the optimal power 
level monotonically decreases to the level ^ and remains at that level afterwards. Similarly, 
the temperature level rises to Tc and remains at that level afterwards. Note that the throughput 
and the energy consumption in Fig. [^ are higher with respect to those in Fig. [51 Parallel to this 
observation, the monotone decrease is sharper in the power policy in Fig. [^ compared to that in 
Fig. [51 Since the power level has to be stabilized at the temperature increase cost paid for 
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Fig. 5. Power, energy and temperature plots for unlimited energy and D = 2 for the single epoch case. 

achieving certain throughput is minimized if energy consumption starts faster and drops later. 

In Fig. |71 we set the deadline to D = 3.5 and the energy limit to E = 17.71. Note that this 
energy level is slightly less than the energy of the power policy in Fig. 0 which translates into 
a right shift of the point to- In particular, we calculate to = 3.2 as the solution of (I5TI)-(I5^ in 
this case. Similar to the effect of decreasing the deadline observed in the comparison of Figs. [5] 
and we observe that decreasing the energy level yields a smoother power policy. Power level 
drops to ^ and the temperature hits at a later time to and both remain constant afterwards. 

In Fig. m we consider the same system as in previous figures with two energy arrivals instead 
of one and with D = 5. In particular, Eq = 6.08 is available initially and Ei = 14.55 arrives at 
time Si = 1.5. In this case, we calculate to = 3.9. The energy causality constraint is tight and the 
power level makes a jump at the energy arrival instant. Note that the temperature is continuous 
at the energy arrival instant even though its first derivative is not. While the power level has a 
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Fig. 6. Power, energy and temperature plots for unlimited energy and D = 3.5 for the single epoch case. 

smooth start, a sharper decrease is observed towards the end since the harvested energy has to 
be fully utilized. In particular, the temperature increase before the energy arrival is kept to a 
minimum level so as to have a higher heat budget for the larger energy that arrives later. The 
temperature hits Tc at f = 3.95 after which the power and temperature both remain constant. 

Finally in Fig. |9l we illustrate a curious behavior in the optimal policy. For this example, we set 
a = 0.1, b = 1.1, Te = 37 and = 37.92. Initial energy is Eq = 25 and energy arrives at f = 2 
with amount Ei = 17 and the deadline h D = 3.5. We observe that energy causality constraint 
is tight at f = 2 whereas it is not tight dXt = D meaning that some energy is wasted in order not 
to cause excessive heat. The temperature generated in this throughput optimal power policy first 
monotonically increases, hits Tc at t = 1.31, remains there till t = 1.66 and drops below Tc. We 
interpret the drop in the temperature in the first epoch as an effort to create temperature margin 
for the high energy arrival in the next epoch. We calculate = 2.23 as the time after which 
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Fig. 7. Power, energy and temperature plots for limited energy E = 17.71 and D = 2 for the single epoch case. 

power level remains at ^ = 10.12 and the temperature remains at Tc. Note that under unlimited 
energy, temperature would hit Tc at t = 0.878. Due to the energy searcity in the first epoch, 
temperature hits Tc later and drops below Tc. A common behavior we observe in each numerical 
example is that temperature ultimately increases between two epochs where energy causality 
constraint is tight. Further research is needed to quantify the relations between the amount of 
temperature generated while performing optimally in terms of throughput. While mono tonicity 
of the temperature is lost when multiple energy harvests exist, we note that monotonicity of the 
temperature is guaranteed in the last epoch due to Lemma [TVl 

VIII. Conclusions 

We considered throughput maximization for an energy harvesting transmitter over an AWGN 
channel under temperature constraints. We used a linear system model for the heat dynamics 
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Fig. 8. Power, energy and temperature plots for two energy arrivals, Eq — 6.08 and Ei — 14.55 al t — 1.5 and D = 5. 

and determined the throughput optimal power seheduling poliey under a maximum temperature 
eonstraint by using a Lagrangian framework and the KKT optimality eonditions. We determined 
for the single energy arrival ease that the optimal power policy is monotone decreasing whereas 
the temperature is monotone increasing and both remain constant after the temperature hits the 
critical level. We then generalized the solution for the case of multiple energy arrivals. While 
monotonicity of the temperature is lost when multiple energy harvests exist, we observed that 
the temperature ultimately increases while maximizing the throughput. We also observed that 
the main impact of the temperature constraints is to facilitate faster energy expenditure subject 
to energy causality constraints. Additionally, even though using all of the available energy is 
optimal for throughput maximization only, with temperature constraints, energy may have to be 
wasted in order not to exceed the critical temperature. 
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Fig. 9. Power, energy and temperature plots for two energy arrivals, i5o = 25 at f = 0 and i5i = 17 at f = 2 and D = 3.5. 
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